On the dangers of using the growth equation on large scales in the Newtonian gauge. 
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We examine the accuracy of the growth equation 5 + 2HS — 4nGpS — 0, which is ubiquitous 
in the cosmological literature, in the context of the Newtonian gauge. By comparing the growth 
predicted by this equation to a numerical solution of the linearized Einstein equations in the ACDM 
scenario, we show that while this equation is a reliable approximation on small scales {k >h Mpc~^), 
it can be disastrously inaccurate (~ 10*%) on larger scales in this gauge. We propose a modified 
version of the growth equation for the Newtonian gauge, which while preserving the simplicity of 
the original equation, provides considerably more accurate results. We examine the implications of 
the failure of the growth equation on a few recent studies, aimed at discriminating general relativity 
from modified gravity, which use this equation as a starting point. We show that while the results 
of these studies are valid on small scales, they are not reliable on large scales or high redshifts, if 
one works in the Newtonian gauge. Finally, we discuss the growth equation in the synchronous 
gauge and show that the corrections to the Poisson equation are exactly equivalent to the difference 
between the overdensities in the synchronous and Newtonian gauges. 



I. INTRODUCTION 

The past few decades have witnessed a remarkable im- 
provement in the level of precision in cosmological obser- 
vations. From measurements of the temperature fluctu- 
ations of the CMB, along with the total energy budget 
and flatness of the universe, to distance measurements 
and determinations of the Hubble parameter, errors on 
the order of a few percent are now commonplace. Within 
an arena of such precision, theorists must be especially 
cautious of deeply embedded approximations within their 
calculations. These approximations may be harmless un- 
der a given set of assumptions (for example, small scales 
and small redshifts), but dangerous if implemented in sit- 
uations where those assumptions are not applicable (e.g. 
future probes on large scales and high redshifts). 

With this in mind, in the present work we examine 
the growth equation in the Newtonian gauge (we discuss 
the relation between the conformal Newtonian gauge and 
the synchronous gauge in Sec|V|, the derivation of which 
involves a number of approximations. The most impor- 
tant of these, as we demonstrate, is the Poisson equa- 
tion which is used to link the matter perturbations (and 
subsequently their growth) to the metric perturbations. 
We show that in linearized general relativity, the Poisson 
equation follows from one of the Einstein constraint equa- 
tions, in which two terms have been discarded. While on 
small scales these terms can be safely neglected, we show 
that on large scales or large redshift, at least one of them 
can be on the same order as the perturbative variables, 
and its absence can introduce significant errors in the 
growth equation. That these terms may be dominant 
on large, horizon size scales has of course been known 
for a long time. What has not been appreciated is that 
on sub-horizon scales these terms may generate contribu- 
tions which are on the order of experimental precision. 
We estimate the error arising from the growth equation 
on various scales and redshifts, and show that the er- 
ror can be significant enough in the context of current 



and proposed experimental limits that caution should be 
exercised in its use in calculations. 

Being aware of this error, and thus avoiding it, is im- 
portant for many models which attempt to break the de- 
generacy between modified gravity (MG) and dynamcial 
dark energy (DDE) (see for example [J 0, i, i, i, [1]). 
These MG studies attempt to discern a deviation from 
general relativity (GR) in observations such as gravita- 
tional lensing, large scale structure growth, and the Inte- 
grated Sachs- Wolfe (ISW) effect. The Newtonian gauge 
is used in these investigations due to the metric pertur- 
bations being the gravitational potential(s) (these two 
perturbations are equivalent in the case of no anisotropic 
stress, and the difference between the two is parameter- 
ized in MG by what is called the gravitational slip). For 
example in lensing studies, it is the gradient of these po- 
tentials which is relevant (which is just given by the Pois- 
son equation in the case of no anisotropic stress), in the 
ISW it is the time variation of the potentials which enter 
the calculation, and in structure growth it is the Poisson 
equation which is used to derive the growth equation. 
Working in the Newtonian gauge then, the Poisson equa- 
tion arises and plays a central role, and is used through- 
out the MG literature. 

In this work we will take a more detailed look at a few 
of these recent investigations of MG, as well as proposed 
measurements at large redshift. As stated above, these 
models, are particularly sensitive to any adjustments to 
the Poisson equation as they introduce modifications to 
the Poisson equation in order to investigate non-standard 
cosmology. The fact that modifications to GR are pro- 
posed to operate at large scales (due to the degree of 
accuracy to which GR has been tested at small scales) 
is also an indication that one should be mindful that 
one is not mistaking a signal for MG when in fact it is 
a sign that one's approximation is becoming less accu- 
rate at the scale of interest. In other words, GR may 
be erroneously thought to break down when it is the ap- 
proximation used, which discards terms inherent in GR, 
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that is the culprit. 

The outhne of our paper is as follows. Except in the 
last section, our entire work is in the Newtonian gauge. 
In Sec|ll] we review the derivation of the growth equa- 
tion and the Poisson equation paying attention to the 
assumptions and approximations that enter the calcu- 
lations. In Sec lIIII we compare the growth predictions 
from the growth equation to a direct numerical integra- 
tion of the linearized Einstein equations in the ACDM 
scenario, and demonstrate that the chief source of error 
in the growth equation arises from the Poisson equation. 
Sec lIVI studies the implications of using the growth equa- 
tion in recent efforts to discriminate general relativity 
from modified gravity theories. Sec|V] deals with relating 
the growth equation and matter perturbations in the syn- 
chronous and Newtonian gauges. Conclusions are found 
in SecEIl 



II. LINEARIZED EINSTEIN GRAVITY: THE 
POISSON AND THE GROWTH EQUATIONS 

In this section we review the derivation of the growth 
equation, highlighting the assumptions and approxima- 
tions that go into the derivation. We work in ordi- 
nary general relativity (GR) with a general stress-energy 
tensor (this can include DDE). We use the notation of 
0,i,i,[3- The perturbed metr ic is of the form 



der, 



-a^ [(1 + 27/^)^, 



2adtBdx'dt 



(1) 



where A^Bjip and E represent metric perturbations and 
a represents the cosmic scale factor. 

Working in Newtonian gauge (we are working with 
coordinate time, and therefore we are not working in 
the "conformal" Newtonian gauge, as we do in SecfVl 
when comparing with the synchronous gauge formula- 
tion), which corresponds to a transformation to a frame 
such that B = E = 0, the gauge-invariant variables char- 
acterizing the metric perturbations become: 



$ = A 

dt 



a^{E + B/a) 



-ilj + a^H{E + B/a) 



A. 



(2) 
(3) 



The energy-momentum tensor can be decomposed as 



-{p + 5p) , 



^ a 



(4) 



where 11^ is a tracefree anisotropic stress, and p, and 
p denote the pressure and density of the cosmic fluid 
respectively. For the moment we take the anisotropic 
stress to be zero. 

The perturbed Einstein equations yield, at linear or- 



+ * = (5) 
-$ + 3i72$ + 3i?$ = -4t:GSp (6) 
H<P + ^ = 4:TrGa{p + p)v (7) 



3$ + 9i7$ 



+ (6H + 6H^ + —)^ = 4ttG(6p + 36p) (8) 
Sp + '6H{6p + 6p) = {p + p) (3$ 

(9) 



[a\p + p)v]' 
aHp + p) 



a 



-(* + ^) (10) 
a \ p + p; 



where a dot, bold or otherwise, denotes a derivative with 
respect to coordinate time t, H = d/a, and any quantity 
preceded by S denotes a perturbation in that quantity. 
These equations determine the two metric perturbations 
$ and along with the velocity potential v and the mat- 
ter perturbation 6p which is conventionally re- written as 
6 = Sp/p. 

We have used the relation between <I> and ^I* in Eq. ([5]) 
in the subsequent equations. This relation simply reflects 
our assumption of no anisotropic stress. From this point 
on we use the Newtonian potential $ to characterize the 
metric perturbation. 

We now switch to Fourier space, an extremely conve- 
nient transformation in the theory of linear perturba- 
tions. In the Fourier-transformed fc-space, all pertur- 
bative variables (e.g. 6p) are replaced by their Fourier 
transforms (e.g. Spk)- For convenience, we suppress the 
k-subscripts and in what follows (unless otherwise men- 
tioned), all perturbative variables can be assumed to be 
in Fourier space. A subscript denotes the present (zero 
redshift) value of that quantity. 

Eq. ([6]) allows us to express the energy overdensity 6 
in terms of the potential $ and the background variables 
as follows: 



inCSp = ^$ + 3i?2$ + 3i7$ 



(11) 



One can see from this equation that on small scales where 
k"^ /a^ ^ the first term on the right hand side dom- 
inates. If one also makes the assumption that the grav- 
itational potential <& is slowly varying (as in the case of 
matter domination) , then one recovers the familiar Pois- 
son equation 



inGSp — $ 



(12) 



These approximations feed into the derivation of an 
equation governing the growth of structure in linearized 
gravity, commonly called the "growth equation" , which 
we now derive. 
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Starting from Eqns.([9]) and ^TO\i we find 



fc2 



-3i7e + 3(l + w)$- (l + w)— w (13) 



-vH{l - 3w) - 



1 + w 



e 



1 + w I + w 



(14) 



where w is the EoS function defined in terms of pressure 
and energy density as w = p/p, O is the combination 
{dp/6p-w)6 [13]. It can be verified that these equations 
are identical to Eq. 30 in [llj . 

The equation which determines the growth of the mat- 
ter perturbation is found to be 



S + S(2- 3w)H + —wS + —(1 + w)$ = 
a 

3(1 + w)[^ + 6(2 - 3w)H] + 3w$ 
-— + 



(15) 



-3i?e 



-(i + 9w)]e 



(one way to see this is by substituting v from Eq. ^ 
and 5 from Eq. pT|) into Eq. (fTU]) . and using the time 
derivative of Eq. ^). 

To get from Eq. [TSl which is exactly true in linear 
order, to the growth equation one must make several as- 
sumptions. First we assume the metric perturbation $ 
is slowly varying ad discard all its derivatives. Next, we 
take 6 ~ 0, 9 = 0, which set the effective sound speed of 
the cosmic fluid dp/ dp approximately equal to its equa- 
tion of state w, and make both of them constant in time. 
We then take and w ~ 0. These approximations are 
perfectly reasonable in a universe heavily dominated by 
pressureless perfect fluid matter, and they lead to the 
following much simplified equation: 



6 + 2H6 



r$ = 



(16) 



Finally using the approximation Eq. (jl2p we recover the 
familiar growth equation: 



S + 2H6 - AnGpS = . 



(17) 



The growth equation is ubiquitous in cosmology and 
is used in many different forms, usually as either a first 
or second order differential equation in a certain growth 
variable, which is a function of d. We list three versions of 
this equation, which will be relevant to this work, below. 

One choice of the growth factor is G = d\n{d / a) / dlna 
which leads to the equation: 



dG 

dlna 



+ 3 



1 dlni?2 

2 dlna 

1 dlni/2 

2 dlna 



G + G^ 



-ilmia) = 



(18) 



A slightly different choice of growth factor is / 
dlnS/dlna leading to the equation 



df 
dlna 



+ f + 



dhiQ„j 
dlna 



f-^GNnra=0 (19) 



Clearly / = G+1. 

A third choice of growth factor 
6„i{k, z)/d„i{k, z — > oo) yields the equation: 



D 



d2 InD 
d (Ina)^ 



1 dH \ dlnD AttG 



H dlna dlna 



pD = (20) 



The question is then: how accurate is the growth equa- 
tion upon use of these approximations? Certainly in 
the matter dominated era {z ^ 1) and on small scales 
(fc ^ .OOlh Mpc~^) these approximations appear justi- 
fied. However, in this era of precision cosmology one must 
be careful to avoid throwing away terms which may ac- 
count for measurable deviations. We now address the 
size and seriousness of such errors. 



III. TESTING AND IMPROVING ON THE 
GROWTH EQUATION 

In this section, we discuss in detail the errors that arise 
in the growth equation from assuming the Poisson equa- 
tion to be true in the Newtonian gauge. We assume a 
model consisting of perfect fluid dark matter and a cos- 
mological constant dark energy. We test the growth equa- 
tion Eq. (|17p by comparing it to a numerical integration 
of the full set of linearized Einstein equations governing 
the growth of perturbations in this ACDM Universe. We 
also ignore radiation in the calculation presented here, 
but we have verified that including radiation does not 
change our conclusions. 

Note that a more complicated model including baryons 
and neutrinos can introduce further sources of divergence 
between the growth equation and the true growth in the 
Newtonian gauge. Baryons have been shown to signif- 
icantly affect the growth of the overdensity (contribut- 
ing an error of 10% on scales below 10 Mpc), and on 
large scales the neutrino anisotropic stress cannot b e ig - 
nored. Details on both these effects can be found in [13 |. 
Stochastic corrections to 5 as a result of random forces 
(possibly arising from the "graininess" of the underlying 
system of particles) was explored in [l^ where it was 
shown that these random forces can lead to significant 
deviations from the nonstochastic solution at late times. 
However, we ignore these effects here, and focus only on 
the errors arising from relativistic corrections to the Pois- 
son equation. 

The zero-th order Einstein and Euler equations de- 
scribing the evolution of the background of this system 
are the following: 



2H + 



SvrGpA 
-3Hp 



(21) 
(22) 



4 



subject to the constraint 



8ttG {p + pa) 



(23) 



Here pA is the energy density of the cosmological con- 
stant. These equations can be solved analytically to 
yield: 



ait) 



1/3 



sinh2/3 ( - 



t 



Hit) = -^cothf-^ 



Pit) 



PAsinh ^ ( ^ 



(24) 
(25) 
(26) 



where tA = 2/ \/2AiKTpA 

The first-order equations governing the evolution of 
perturbations can be easily derived. Using vj = —va we 
find from equations dH), ^ and pi)) : 



-4F$ - SttGpa* 



Equations 



5 = 3$+ 

Vf = -$ 

and ([7]) yield the constraints: 

/ \ P 

3i7 + +— $ = -47rG(5p 



(27) 

(28) 
(29) 

(30) 
(31) 



An alternative method of deriving Egs.^ ipTjl - ipTI) ) is to 
start with the zero-th and first order Einstein equations 
for a Universe consisting of matter and scalar field dark 
energy (these can be found in, for instance [Tsj in the 
Newtonian gauge and in the synchronous gauge), 
and then to take the limit as the scalar field tends to 
a cosmological constant, that is, the potential V {(f)) — 
PA — constant and (t){t — Q) = Q. 

For the numerical integration, we set initial conditions 
at a redshift z — 1100 where we set 5 = 10~^ and $ = 0. 
Equations Eq. ([30| and Ea. ([3T)) are used as independent 
checks on the accuracy of our numerics. We call the 
growth predicted by our numerical evolution the "true 
growth" . A simple mathematica notebook which per- 
forms the above numerical integration can be found at 
0. 

Henceforth, we denote the true growth by the variable 
5 and that predicted by the growth equation pT|) by 5g. 
We express the percentage departure of 5g from 5 by the 
quantity 



(^9 - 5) 



(32) 



Fig. [T] shows A as a function of scale, and we can 
see that the growth equation fails spectacularly as one 
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FIG. 1: The (percentage) error A in the growth ((5gi) pre- 
dicted by the usual growth equation Eq. (|17|) as a function of 
scale, for three different redshifts 
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FIG. 2: The (percentage) error A in the growth [Sg) predicted 
by the usual growth equation Eq. (|17|l as a function of redshift, 
for four different scales. From the top down, the curves are 
for k = .OOlh Mpc-\ k = .Olh Mpc-\ k = .Ih Mpc-^ and 
k —h Mpc~^ respectively. 



approaches the horizon scale. This is also reflected in Fig. 
[21 which shows A as a function of redshift for different 
scales. The failure at horizon scales is not a surprise of 
course. What is of interest is that there are significant 
deviations arising at scales within reach of experiment, 
and at levels that could be misinterpreted as a signal 
of the breakdown of GR if one were unaware that the 
Poisson equation approximation itself is breaking down 
at these scales. 

Aside from ACDM, we have also verified the failure of 
the growth equation in a model where the dark energy is 
a scalar field with a quadratic potential (with mass com- 
parable to the present value of the Hubble parameter), 
e.g. the scenario considered in detail in p^ . 

Having established the failure of the growth equation 
on large scales, we now proceed to investigate the cause 
for this failure and propose an improved version of this 
equation. On large scales, the most obviously suspect 
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step in the derivation of Eq. (jT?!) is the replacement of 
Eq. (|lip by the Poisson equation, or in other words, the 
neglect of the term 3H{H^ + $) in comparison to the 
term (^k^/a^) $. While this is justified on small scales 

^ a?H'^, it is not justified on scales close to the hori- 
zon, where the terms 37f^$ and [k'^/a'^) $ are on the 
same order. The error introduced in neglecting a term 
of the same perturbative order is rapidly magnified in 
the process of integrating over the age of the Universe, 
leading to the gigantic deviation from the true growth as 
shown in Figs[T][21 

A smaller error is introduced in ignoring the velocity 
of the metric perturbation $. However, we have verified 
that this does not lead to a significant error except on 
horizon scales. 

To understand the magnitude of the error introduced 
by ignoring the term proportional to let ^ (z, k) denote 
the ratio of the 3iJ^$ term to the (fc^/a^) $ term, i.e.. 



(33) 



We can estimate the size of this term (in the matter dom- 
inated regime, for simplicity) by using 



l + z 



(34) 



where the subscript zero denotes today's value, and oq = 
1. For convenience, we express k as 



k = 10""- 



Mpc 



(35) 



such that n now denotes the scale. Finally, the present 
value of the Hubble parameter Hq is expressed the usual 
way, as 



Ho = 100 



h km 
Mpc s 



(36) 



Plugging the above into Eq. ([55)) . using r2,„o — -3 (and of 
course restoring the correct units by dividing by the speed 
of light squared), we obtain the following simple formula 
for ^ which shows its dependence on both physical scale 
and redshift: 



e= (l + z)102 



(37) 



In Figl3] we plot ^ as a function of z at various scales. 
One can see that the discarded term {3H^^) can be a 
substantial, non-negligible fraction of the retained term 
(^k^ /a^) $ as n becomes larger than 2 (i.e., k < lO^'^h 
Mpc~^), or at large redshifts. At n « 3, even at small 
redshifts this term will be > 0(10)%. As mentioned be- 
fore, the cumulative effect of integrating over the history 
of the Universe causes the error to magnify rapidly. 

Based on the above considerations, we propose the fol- 
lowing modification to the growth equation 



1 + ^ 



(38) 
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FIG. 3: The correction term ^ as a function of redshift, for 
different scales according to Ea. (|37|l as function of redshift, 
for four different scales. From top down, the curves are for 
k = .OOlh Mpc"\ k = .Olh Mpc"\ k = .Ih Mpc"^ and fc = ft 
Mpc~^ respectively. 



Figs|3][5l which show the A arising from the improved 
growth equation (let us call this 6gi), verify that Eq. ([551) 
provides a far better approximation to the true growth 
than the usual growth equation. Even on scales of the 
size of .Olh Mpc~^ , the error is on the order of a percent 
for all redshifts. For scales on the order of the horizon, 
the error is somewhat large for low redshifts, presumably 
as a result of the dark energy dominance which causes 
the $ term in Eq. [Tl]to become significant. 

The cumulative effect of integrating a small error over 
a large time period can be demonstrated by considering 
the difference between the usual growth equation Eq. p7p 
and the improved growth equation Ea. (j38p . Taking the 
growth predicted by the improved growth equation Sgi 
as a proxy for the true growth (according to FiglH the 
difference between S and Sgi is negligible up to scales of 
.Olh Mpc"^), one can define the error variable in terms 
of Sgi (rather than S as in Eq. ((32)) 1: 



(39) 



where 6g and Sgi evolve according to Eqs. lfTT]) and ((SSj) 
respectively. The growth variables defined in this section 
are listed for convenience in Table HI 

One can now construct a differential equation for the 
evolution of A for a given k (under the assumption of 
matter domination): 



1+e 



1 + e 



1 + A 



(40) 



Evolving this equation with the initial conditions 

A{ti) — A{ti) — 0, one obtains the almost identical re- 
sults as shown in Fig. [2] for scales up to .Olh Mpc~^. In 
other words, this shows that the large error in Figs. [1] 
and [5] are (largely) the result of ignoring the 3-ff^$ term 
in comparison to the (fc^/a^) $ term for large scales. 
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Variable 


Definition 


S 


growth according to numerical integration of equations (I27I)-(I29I) 




growth according to the usual growth equation Eq. 1171 


Sgi 


growth according to the improved growth equation Eq. |38l 


A 


{Sg - 5) /S in Figsd] and [land {Sgi - 5) /S in Figs|4]and[5] 


A 


{Sg - Sgi) /Sgi 


G 


dln{S / a) / dlna 


f 


dlnS 1 dlna 



TABLE I: List of growth variables defined in the text. 
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FIG. 4: The (percentage) error A in the growth (Sgi) pre- 
dicted by the improved growth equation from Eq. (I38p as a 
function of scale, for three different redshifts 



We wish to emphasize that contrary to what is often 
suggested in the literature (e.g. 0), the growth factor 
is not scale independent. This can clearly be seen from 
the form of the ^ correction Eq. ((33|) . Hence an observed 
scale dependence in the growth factor cannot be taken to 
imply a departure from Einstein gravity. 

Finally, note that the other versions of the growth 
equation mentioned previously need to modified by such 
a term as well. For example, in terms of the growth factor 
G, the growth equation Eq. p^ becomes 



dG 
dlna 



1 dlni?2 
4-1 

2 dlna 



^ 2 dlna ~ 2 1 + ^ 



(41) 



Based on the above discussion, we recommend the use 
of Eq. ((38| over the usual Eq. (|17p as a much better ap- 
proximation to the linear growth of matter perturbations 
on large scales. The importance of using a version of the 
growth equation accurate to large scales and high red- 
shifts is underscored by the fact that scales on the order 
of .Olh Mpc~l and redshifts > 1 are within the reach 
of future surveys such as ADEPT [l^ as well as surveys 
based on the 21 centimeter emission. 



FIG. 5: The (percentage) error A in the growth iSgi) pre- 
dicted by the improved growth equation Eq. (I38|l as a func- 
tion of redshift, for four different scales. From the bottom up, 
the curves are for k = .OOlh Mpc"\ k = .Olh Mpc"\ k = .Ih 
Mpc~^ and k = h Mpc~^ respectively. 



IV. PROBING BEYOND-EINSTEIN PHYSICS 

We now examine the implications of the failure of the 
usual growth equation on some recent studies which have 
used the growth equation as a probe to distinguish Ein- 
stein gravity from new physics. We first consider the 
work by Linder |3i| and Linder and Cahn In 
the author finds that the linear growth of structure ac- 
cording to Einstein gravity (represented by Eq. pjp ) 
can be modeled with remarkable accuracy by a single 
parameter 7 ~ 0.55 for the entire expansion history 
of the Universe. In particular, if one uses Eq. (flS)) . 
then the growth history, characterized by the quantity 
G(z) = d\n{S/a) /d\na can be described by the rela- 
tionship 



G(a) = flm (a) 



1 



(42) 



The parameter 7, dubbed the "growth index", is given 
by the fitting formulas 



7 
7 



0.55 
0.55- 



.02 [1 
.05 [1 



w [z 
w [z 



1)], 
1)], 



w < 
w > 



(43) 
(44) 



over the whole range < a < 1. 



7 



In [3| the authors provide analytical arguments sup- 
porting the above formulas, and contend that the narrow 
range in which this parameter is constrained in the con- 
text of Einstein gravity, allows for a possible distinction 
of Einstein gravity from other gravitational theories in 
which this parameter is not constrained in this range. It 
is not difficult to repeat the steps of their analytic calcu- 
lation, but now with the ^ term included. One can solve 
Eq.(|4T]) formally 



G(a) = -1 + 



(45) 



where the substitution Q,„ = 1 ~ ^w, where fl^ is 
the dark energy fraction, has been made. In the mat- 
ter dominated regime we use the fact that flm S> ^Iw 
The term can be neglected, and using the form for 
the Hubble parameter in a matter dominated universe: 



a •'(1 -I- ri^(a)/f2m(a)), G is found to be 



G(a)^-I + ^(l-^ 



(46) 







a' 



a' 



2(1 + 
3 



2 4(1 + 



The corrected, ^ dependent, growth index can then be 
recovered by using 
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G{a) 
fitu(a) 



(47) 



in Eq.dUD) 

To see how well the fit described in equations 
(|44)) works for the true growth (as opposed to the growth 
predicted by the Eq. ([TT])) we plot the growth function 
/ = d In In a = G + 1 as derived from the true growth 
(on the plots this is denoted by f{z;d) ) against the fit 
(a)°'^^ for the redshift range < z < 10. The growth 
function derived from the growth equation (denoted by 
^ is shown for reference. 
Fig. [HI confirms the findings of [1, |3| for the growth 
equation (and hence for small scales) , but shows a drastic 
departure of the 7-fits from the true growth for scales 
k < .Olh. For these large scales, it turns out that 7 has a 
very strong z-dependence, and the simple fits of Eq. 
(|44p do not refiect the true growth, even approximately. 

For the true growth on large scales, behavior of the 
growth index 7 = ln(/) / ln(f2m) is shown for three differ- 
ent scales k = Ah Mpc~\fc = .Olh Mpc^^ and k = .OOlh 
Mpc~^ in Fig. [71 It is clear that for the true growth on 
large scales, 7 has a strongly non-linear z dependence and 
is not by any means restricted to a small range. Hence 
on these scales it cannot be used to discriminate between 
GR and modified gravity theories. 
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FIG. 6: The growth parameter /(z; S) derived from the true 
growth, plotted against the numerical fit fim [z)'^^ for three 
different scales. The growth parameter j(z\ Sg) derived from 
the growth equation is shown for reference. 



Polarski and Gannouji [s"] have proposed an interest- 
ing strategy for discriminating GR from modified gravity, 
using the present values of 70 = 7(2; = 0) and its deriva- 
tive 7q ^ j'{z = 0). Starting with the growth equation in 
the form of Eq. (jl9p . they obtain the following constraint 
condition finking 70, 7o, O^.o, and WeS,o = wde,o^de,o 
(where wde is the equation of state parameter of the dark 
energy component): 
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FIG. 7: The growth index 7(2:) derived from the true growth, 
plotted against redshift for three different scales. Contrary 
to the analysis in [3, 01 , 7 seems to have a strong and highly 
non-linear z-dependence on large scales. 



7o 



Inn 



m,OJ 



-0.''° 



3(70 - ^)Woff,0 



m,0 



(48) 



Based on this equation, Polarski et. al. then derive the 
constraint |7o| < 0.02 for ACDM. They therefore con- 
clude that a measurement of 7g outside this range would 
signal a departure from ACDM. 

However, in the hght of our discussion in Sec lIIIi the 
correct starting point should be the modified growth 
equation Eq. ([55]) . This leads to the following modified 
form of the constraint equation: 



7o = (Inf^mlo) 



-^1% - 3(70 - ^) Wo 



O "m,0 
'2 1+f 



(49) 



Clearly, the ^-correction can cause the value of 79 
to depart from the range observed by Polarski et. al. 
quite drastically. As an example, consider the param- 
eter choices wde = —1, ^m.o = 0.3 and 70 = 0.554. 
Eq. yields 7q = —.0192, which is within the range 
mentioned above. For scales of fc > O.lh Mpc""'^ and 
smaller, Eq. agrees. However, for k — O.Olh Mpc~^, 
and k = O.OOl/i Mpc-\ Eq. dH]) yields 7^ = -0.040 
and 7o = —0.565 respectively, indicating that the bounds 
derived using the growth equation are not respected on 
large scales. 

Acquaviva et. al Q have suggested a null test param- 
eter e (k, a) = flm (a)"'''''^'' / (a) — 1 (where j{z) is given 
by Eq. (|44p) as a tool to discriminate GR from modified 
gravity . They claim that any non-zero measurement of e 



FIG. 8: The null test parameter e(fc, 2) as proposed in 0] as 
a function of redshift for different scales. 



which cannot be attributed to systematics should be in- 
terpreted as an signature of modified gravity. However, 
it is clear that the e parameter is based on the growth 
equation Eq. pT|) . and given the failure of the growth 
equation on large scales, this parameter is not reliable. 
To demonstrate this point, we plot the value of this pa- 
rameter vs redshift in Fig. [5] for three scales. We see that 
on small scales {k ~ O.lh Mpc~^) the parameter is close 
to zero (as expected) , but strongly deviates from zero for 
large scales and large redshifts. 



V. CHOICE OF GAUGE 

In this section examine the relationship between the 
growth equation in the synchronous and conformal New- 
tonian gauges (so far we have not used the conformal 
time form of the line element, but we will here for eas- 
ier comparison with previous works) . We rely heavily on 
the work of Ma and Bertschinger , and any notational 
differences between our work and theirs will be specified 
for clarity. The line element in the synchronous gauge is 
given by 



a'^{T){-dT'^ + (5y + hij)dx^dx^) 



(50) 



The Einstein equations are written in terms of two func- 
tions, h{k, t) (which is the trace of hij) and ri{k, r), where 
k = kk, that are defined through the Fourier integral of 
the metric perturbation hij as 

h,j{x,T) = I d^ke'''-'^[kikjHk,T) (51) 



+ ihkj - -Sij)6r]{k,T)] 

In pr) ] the line clement in conformal Newtonian gauge is 
written as 

ds^ = a^{T)[-{l + 2^P)dT^ + (1 - 2(l})dx'dx,] (52) 
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In the case we are examining where anisotropic stress is 
absent we can relate the variables (/) and ip to our $ as 

= -0 = 'J'- 

The variables in the synchronous and conformal New- 
tonian gauges can be shown to be related by 

= -^{h"ik,T)+6ri"{k,r) (53) 
a 

where the prime denotes a derivative with respect to 
conformal time r, whereas a dot represents a derivative 
with respect to coordinate time t. Using the variable 
a = {h' + 6ri')/2k'^ this relation can be written as 

$ = ad + da = ^(aa) (54) 

The Einstein equations in the synchronous gauge in 
the absence of anisotropic stress are given by 

k^n-l-h' ^ -AnGNa'^Sp, (55) 
2 a 

k^T]' = 4:TrGNa^{p + p)e (56) 

h" + 2—h' -2k^ri = -SnGNa'^iiSps) (57) 
a 

h" + 6ry" + 2-{h' + Qr^') - 2k^T] ^ (58) 
a 

where 9 is related to the fluid velocity by a divergence, 
9 = ik^Vj, and the subscript s denotes the synchronous 
gauge (we follow [ill here and use c to denote the con- 
formal Newtonian gauge as well). In the simple case of 
a matter dominated universe one can derive the growth 
equation for = 6ps/p from these equations, and one 
arrives at 

6s + 2H5s - At^Gn5s = (59) 

We see that the growth equation is exact in the syn- 
chronous gauge in the case of matter domination. In 
order to compare this to the case of the conformal New- 
tonian gauge we make use of the relation between the 
perturbations in the two gauges 

5c- a- ^ 5c + '3Haa{l + w) (60) 
P 

(In the matter dominated regime one can obviously sim- 
plify this to ds = Sc + 3Haa). Inserting this relation into 
the growth equation in synchronous gauge, Eq. (j59p . and 
using the relation Eq. (|54p we flnd 

S'c + 2Hdc - 'inGNSc - 3-ff^$ + 3ff$ = (61) 

Comparing this to the relation we found for 5 in Eq. (jl5p 
we see that they exactly match as one would expect. In 
order to see this starting from Eq. pS]) , one should use the 
approximations of matter domination and set w, ?«, 0, 
to zero, and use Eq.® to substitute for $. 



What this shows is that if one wishes to use the exact 
growth equation, one should work in synchronous gauge. 
The difference between the growth equation in the two 
gauges is given by the extra $ and $ terms as we have 
noted, which cause a large deviation in the conformal 
Newtonian gauge from the solution to the usual growth 
equation. In other words, the deviation from the usual 
growth equation in conformal Newtonian gauge is a mark 
of the deviation between the evolution of the matter per- 
turbation in the synchronous and conformal Newtonian 
gauges. 

In the process of observation, one is essentially viewing 
a weighted two-dimensional projection d (n) of a three- 
dimensional field S{n) . 6 (n) and 6 (n) are linked through 
a line-of-sight-integral 

POO 

S(n)= / dzWx (z) S {fir (z) , z) (62) 

"'0 

where n is a direction in on the sky, r{z) is the comoving 
distance to a point at redshift z, and Wx (z) is a window 
function which selects a range along the radial coordinate 
that contributes to the observable 6. (For more details, 
see e.g. [H) 

6 (h) is obviously a gauge-dependent quantity, but 
as we have demonstrated above, for subhorizon scales 
smaller than (10~^h Mpc""'^), the gauge choice does not 
lead to a significant difference in 6. It is only on scales 
close to the horizon that the two gauges can diverge sig- 
nificantly. In some recent work [20|, it has been shown 
that the Poisson equation is valid on all scales if one 
works in a comoving-orthogonal gauge. In our opinion, 
the question of which gauge works best in modeling real 
astronomical measurements on large scales is important 
and merits further research. 



VI. CONCLUSIONS 

With the current and future ability of cosmological 
probes to collect data with high and ever increasing 
precision, one must be careful when incorporating ap- 
proximations into calculations which will be compared 
to such data. To this end, we have tested the familiar 
growth equation against direct numerical simulations of 
the growth of perturbations in a ACDM Universe and 
found it to be strikingly inaccurate on large (fc < O.lh 
Mpc~^) scales in the Newtonian gauge. We have traced 
the source of the inaccuracy to general relativistic correc- 
tions to the Poisson equation which become important at 
large scales and at large redshifts within the Newtonian 
gauge. We proposes a modified version of the growth 
equation for use in the Newtonian gauge, which we show 
to be highly accurate on all scales up to a tenth of the 
horizon, for all redshifts. We have examined the implica- 
tions of the failure of the growth equation on recent ef- 
forts to distinguish GR from modified gravity, finally we 
have discussed the growth equation in the synchronous 
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gauge and demonstrated that the corrections to the Pois- 
son equation correspond exactly to the difference between 
the overdensities in the two gauges. Our results have im- 
portant implications on efforts to design null-test param- 
eters to distinguish between models of gravity, as well as 
on constraining cosmological parameters based on obser- 
vations from future probes. As modifications to general 
relativity may begin to operate on large scales, one must 
take care to ensure that approximations are not mas- 
querading as beyond Einstein physics. 
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